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Abstract. This paper studies the one-loop expansion of the amplitudes of electromag- 



netism about flat Euclidean backgrounds bounded by a 3-sphere, recently considered in 
perturbative quantum cosmology, by using (-function regularization. For a specific choice 
of gauge-averaging functional, the contributions to the full £(0) value owed to physical de- 
grees of freedom, decoupled gauge mode, coupled gauge modes and Faddeev-Popov ghost 
field are derived in detail, and alternative choices for such a functional are also studied. 
This analysis enables one to get a better understanding of different quantization techniques 
for gauge fields and gravitation in the presence of boundaries. 
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Gauge-averaging junctionals for Euclidean Maxwell theory ... 
1. Introduction 

The way in which quantum fields respond to the presence of boundaries is responsible for 
many interesting physical effects (e.g. the Casimir effect), and plays a very important 
role in quantum gravity and quantum cosmology. In that case, the (formal) quantization 
of gauge fields and gravitation via Wick-rotated Feynman path integrals is expressed in 
terms of quantum amplitudes of going from a 3-metric and a field configuration on an 
initial spacelike surface to a 3-metric and a field configuration on a final spacelike surface. 
Whilst mathematics enables one to understand which compact boundary geometries do 
actually exist, the methods of quantum field theory fix the boundary conditions for scalar, 
fermionic and gauge fields, as well as gravitation and corresponding ghost fields for spins 
1, | and 2. Although the full theory via path integrals is in general ill-defined, since there is 
little understanding of the measure for quantum gravity and of the corresponding sum over 
all Riemannian 4-geometries, the one-loop approximations of these ill-defined functional 
integrals can be evaluated in terms of well-defined mathematical concepts [1-9] . Since one 
has then to study determinants of second-order, self-adjoint, elliptic operators, the basic 
tool used by theoretical physicists is the generalized Riemann ^-function formed by the 
eigenvalues of these operators as [9] 

oo oo 

= X X d m(n)Xn,m~ S ■ (1-1) 
n=no m=m 

With our notation, n, m are degeneracy labels, and d m {n) is the degeneracy of the eigenval- 
ues A n>m , which is taken to depend only upon the integer n, as happens in many physically 
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relevant applications (including the ones described in this paper). The regularized C(0) 
value yields both the scaling of the one-loop prefactor and the one-loop divergences of 
physical theories. 

In particular, the problem of the one-loop finiteness of (extended) supergravity theories 
in the presence of boundaries is still receiving careful consideration in the current literature 
[1-9]. As emphasized in [9-11], one can perform one-loop calculations paying attention 
to: (1) S-matrix elements; (2) topological invariants; (3) presence of boundaries. For 
example, in the case of pure gravity with vanishing cosmological constant, A = 0, it is 
known that one-loop on-shell S-matrix elements are finite. This property is also shared by 
A = 1 supergravity when A = 0, and in that theory two-loop on-shell finiteness also holds. 
However, when A ^ both pure gravity and N = 1 supergravity are no longer one-loop 
finite in the sense (1) and (2), because the non-vanishing on-shell one-loop counterterm is 
given by [10] 

In (1.2) ?= n — 4 is the dimensional-regularization parameter, \ is the Euler number, S is 
the classical action on-shell, and one finds [9,10]: A = B = — for pure gravity, and 
A= ||, S = — j| for A = 1 supergravity. Thus, B ^ is responsible for lack of S-matrix 
one-loop finiteness, and A ^ does not yield topological one-loop finiteness. 

In the presence of boundaries, however, a much larger number of counterterms can be 
obtained even just at one-loop order in perturbation theory, using the extrinsic-curvature 
tensor and the Ricci tensor of the boundary. Thus, if any theory of quantum gravity 
can be studied from the perturbative point of view, boundary effects play a key role in 
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Gauge-averaging junctionals for Euclidean Maxwell theory ... 
understanding whether it has interesting and useful finiteness properties. It is therefore 
necessary to analyze in detail the structure of the one-loop boundary counterterms for 
fields of various spins, and the techniques developed so far are described in detail in [2-9]. 
The corresponding problems are as follows. 

(i) Choice of locally super symmetric boundary conditions [1-4,8-9]. They involve the 
normal to the boundary and the field for spin |, the normal to the boundary and the spin-| 
potential for gravitinos, Dirichlet conditions for real scalar fields, magnetic or electric field 
for electromagnetism, mixed boundary conditions for the 4-metric of the gravitational field 
(and in particular Dirichlet conditions on the perturbed 3-metric). 

(ii) Quantization techniques. One-loop amplitudes can be evaluated by first reducing 
the classical theory to the physical degrees of freedom by choice of gauge and then quan- 
tizing, or by using the gauge-averaging method of Faddeev and Popov, or by applying the 
extended-phase-space Hamiltonian path integral of Batalin, Fradkin and Vilkovisky [2-9]. 

(iii) Regularization techniques. The generalized Riemann ^-function [12] and its reg- 
ularized £(0) value can be obtained by studying the eigenvalue equations obeyed by per- 
turbative modes, once the corresponding degeneracies are known, or by using geometrical 
formulae for one-loop counterterms which generalize well-known results for scalar fields, 
but make no use of mode-by-mode eigenvalue conditions and degeneracies [2-9,13]. 

Since the various quantization and regularization techniques mentioned so far have 
been found to give rise to different estimates of the ((0) value for all spins > 0, it is crucial 
to get a better understanding of the C(0) values obtained by using the manifestly gauge- 
invariant quantization techniques previously listed. The aim of this paper is to perform 
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this analysis in the simplest (but highly non-trivial) case, i.e. the one-loop amplitudes of 
vacuum Maxwell theory when a 3-sphere boundary is present. One is thus led to make a 
3 + 1 split of the 4-vector potential, expanding its components on a family of 3-spheres 
centred on the origin as [9,14] 

oo 

A (x,r) = J2 R n(r)Q (n) (x) (1.3a) 

n=l 

A k (x,r) = Al(x,r) + A%(x,t) for all fc = 1,2,3 (1.36) 

where Q^ n \x) are scalar harmonics on the 3-sphere, whereas the transverse part A\ and 
the longitudinal part A% are expanded, Wk = 1,2,3, as 

oo 

A T k (x,T) = J2fn(r)Sl n \x) (1.4) 

n=2 

oo 

A L k { X ,T) = Y J 9n{T)Pi n \x) . (1.5) 

n=2 

Of course, the (x) and ( x ) are the transverse and longitudinal vector harmonics on 
S 3 respectively, and their properties are described in detail in the appendix of [14]. Note 
that, strictly, normal and tangential components of A^ are only well-defined at the 3-sphere 
boundary, where r = a, since a unit normal vector field inside matching the normal to S 3 
at the boundary is ill-defined at the origin. As in all mode-by-mode calculations, we are 
performing a local analysis, where one takes that scalar field whose expansion on a family 
of 3-spheres centred on the origin matches the A (x,a) value at the boundary. Moreover, 
one takes that 3-vector field whose expansion on a family of 3-spheres centred on the origin 
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matches the Ak(x, a) value on S 3 . One then has to show that an unique regular solution 
of the corresponding boundary-value problem exists, although the unit normal vector field 
inside is ill-defined at the origin. It will be shown in section 3 that this is indeed the case. 
[We are grateful to Dr. A. Kamenshchik for correspondence about this problem] 

Our paper is thus organized as follows. Section 2 derives the contribution of the 
physical degrees of freedom (i.e. the modes f n (r) appearing in (1.4)) to the C(0) value, 
following [9,14]. Section 3 studies the coupled set of second-order ordinary differential 
equations expressing the eigenvalue equations obeyed by the gauge modes g n and R n , 
Wn > 2. Section 4 derives the contribution of the -Ri-mode of (1.3a), which remains 
decoupled from R n (r) and g n (T), Vn > 2. Section 5 studies the corresponding form of the 
ghost operator, and various possible choices of the gauge-averaging term in the Faddeev- 
Popov formula. Open problems and concluding remarks are presented in section 6. 



2. Physical degrees of freedom 



Within the Faddeev-Popov approach to Euclidean Maxwell theory one deals with gauge- 
invariant amplitudes of the form [15] 



Z[g) = J MA]6mA)] det 



6$(A) 
8e 



exp 



/ -F^F^^Tgd A x 
Jm 4 



(2.1) 



where A M is the 4-vector potential, F^ v = d^A v — d v A [1 denotes the electromagnetic- 
field tensor, and g is the background 4-metric. These amplitudes are more conveniently 
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re-expressed as 



Z[g}= / nM]v2[c, c*] exp ( - Ie) 



(2.2) 



where the total Euclidean action Ie = Ie + Iga + I g h = Ie + I g h is given by 




d>(A) 



a/ det g d A x 



Ie = Igh + 



2a 



(2.3) 



In these formulae is an arbitrary gauge-averaging functional which depends on the 

U(l) potential A and its covariant derivatives, and a is a positive dimensionless parameter. 



a suitable measure on the space of connections, a suitable measure for ghosts, and the 
Faddeev-Popov determinant respectively. The inclusion of gauge-averaging functionals 
and corresponding ghost fields (cf section 5) is necessary to extract the volume of the 
gauge group. In other words, on integrating over all field configurations one integrates 
infinitely many times over the volume of the gauge group, whereas we need to concentrate 
the measure over a subset of configurations containing a single point for each orbit of the 
gauge group. This is achieved using (2.1)-(2.3). 

In recent years, the case of flat Euclidean backgrounds bounded by a 3-sphere has been 
studied as the first step of a program aiming to get a better understanding of one-loop 
properties of supersymmetric field theories in the presence of boundaries [9] . As described 
in section 1, for this purpose one is then led to make a 3 + 1 split of the 4- vector potential 
as in (1.3)-(1.5). In this section we compute the contribution of the physical degrees of 



Igh is the corresponding ghost-field action. Moreover, [ii\c, c*], det 



are 
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Gauge-averaging junctionals for Euclidean Maxwell theory ... 
freedom / n (r) of (1.4) to the one-loop amplitudes of vacuum Maxwell theory in four- 
dimensions. If the measure in the path integral is scale-invariant (see [9,14] and references 
therein) such Z^> amplitudes take the asymptotic form 

Z (1) (a) ~ W a m e' 1 (2.4) 

where W is an arbitrary constant and a is the 3-sphere radius. 

The physical modes f n {i~) are always decoupled from the gauge modes by virtue of 
the properties of the transverse vector harmonics. The corresponding elliptic operator in 
the Euclidean action is found to be [9,14] 

D n = ---^(r4-) + 1 ^ for all n> 2 (2.5) 

T (IT V QjT / T 

whose eigenfunctions are f n (r) = A n J n (\^T), A n being a constant [16]. Since locally 
supersymmetric boundary conditions require that either the magnetic field or the electric 
field should vanish on S" 3 [9,14], one has to compute the regularized C(0) value for the 
generalized zeta- function obtained from the eigenvalues of D n when f n (r) is subject to 
Dirichlet conditions on S" 3 (i.e. magnetic case) or Neumann conditions on S 3 (i.e. electric 
case). In the magnetic case, following the detailed analysis of [14], one finds 

d P£,F, (0) = (2.6) 

where the label (PDF) reminds us that (2.6) is the contribution of the physical degrees 
of freedom to the full C(0). In the electric case, following section 5.9 of [9], we begin 
by taking the Laplace transform of the heat equation, where J n (\^Er) is replaced by 

8 



Gauge-averaging junctionals for Euclidean Maxwell theory ... 
the linear combination A n I n (ar) + B n K n (ar). The ratio ^ is found by requiring that 

-^:^A n I n {ar) + B n K n (ar)^(a) = 0, Vn > 2, which takes into account the eigenvalue 

condition J n (\^Ea) = 0, where a is the 3-sphere radius. Thus, the Laplace transform of the 
kernel of the heat equation for spin 1 when f n (a) = 0, Vn > 2, is an infinite sum of products 

G n of functions G n of the type G n = T (^A n I n (aT) + B n K n (aT)^ . More precisely, defining 
t< = mm(r, r'), t> = max(r,r'), one finds that G n (t, t', a 2 j = G n ^r<, a 2 ^jG n ^r>, c 



where 



G n r<,cr = r</ n (crr<) 



G n T>,(7 = T> 



(2.7a) 



(2.76) 



This implies that the free part of the heat kernel is equal to the one found in [14], and 
hence does not contribute to C(0)- We therefore study the interacting part [9] 



G int I a 2 



n=2 



" " 1 1 in — \tt ( — T/( n <™) 



(2- 



where /(n; era) is the function defined in equation (4.1.18) of [9]. Thus, we have to work out 
the uniform asymptotic expansions of the various terms on the r.h.s. of (2.8) according to 
the relations (4.4.13)-(4.4.22) of [9]. Setting a = 1 for simplicity, and defining y = -^==, 
this yields 



Kn(o-) I'M 



My) + 



Ai(y) , A 2 (y) A 3 (y) 



n 



+ 



+ 



+ ... 



(2.9) 
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where [9] 



A (y) = -1 



A 1 (y) = -y(l-y 2 ) 

A2{y) = - y ^{i-y 2 ) 



(2.10) 
(2.11) 

(2.12) 



My) = -fg (i - y 2 ) (io - esy 2 + 7 V 



(2.13) 



Thus, since f(n; a) ~ Vn2 + (J " 



, g2(y) , g 3 (y) , g 4 (y) , 
n "r IP 1 n 3 1 n 3 r "' 



, using the explicit forms 



of Bi(y) appearing in equation (4.4.12) of [9] one finds 



l^W)^ *~ 



CM C 2 {y) CM CM 



n 



rr 



(2.14) 



where [9] 



d(y) = - y -(l - y 2 ) 
c 2 {y) = V ^(i-y 2 )(iy 2 -i 



C 3 (y) = -^(i-y 2 ) ( 3 -20^ + 21/) 



(2.15) 



(2.16) 



(2.17) 



CM = -^{l-y 2 ){9- 122y 2 + 301y 4 - 196y 6 



(2.18) 



Note that additional terms in (2.14) have not been computed since they give a contribution 
equal to 0(\/i) (t being a parameter not related to Lorentzian time), and hence do not 
affect the £(0) value. Thus, using (2.8) and (2.14)-(2.18) and taking the inverse Laplace 
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transform, one finds that the integrated heat kernel has an asymptotic expansion as t — > + 
given by [9] 

oo 4 

G int (t) ~ - (" 2 - l ) E ^ *) + °(^) ( 2 - 19 ) 



n=2 i=l 



where [9] 

7i(M) = -ie-» at (2-20) 

7 2 („,t) = l-^ = n 2 e-» a «-JIe-» a « (2.21) 

/ 3 (n, = -|te-» a * + ^n 2 e-» a * - Y^e~ nH (2.22) 

/ 4 (n,t) = -^^ 2t + -4^- n2t -— 4« 4 e- at + — 4- 6 ^ • (2-23) 
J V ; 4v/i 15v^F 105 v^F 945 ^ v ; 

The interacting part G tnt (t) is an even function of n, and we can compute its contribution 
to C(0) using the Watson transform defined in [9]. The poles of the integrand at and 
at ±1 are excluded, since the sum over all n in (2.19) only starts from n = 2. The poles 
at ±1 do not contribute, because the integrand of the Watson transform has zeros at ±1. 
In our case the constant contribution arising from the poles is given by the constant term 
appearing in the inverse Laplace transform of — ^/(0,cr), which is equal to |. One thus 
finds 

4 



G tnt (t) ~ 



- + o(Vi)- (>-i)5>(p,t)dp 



7T 3 1 , 55 ,— l 1 1 1 . f~ . 
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where the contributions — | and — are owed to (2.21) and (2.23) respectively. One thus 
obtains the PDF value 

<if DF, (0) = ^ . (2.25) 



3. Coupled gauge modes 



The Euclidean action is obtained multiplying the Lorentzian action by — i, setting t = —ir, 
and bearing in mind that (Aq)l^ = (A^Edr, so that the r n -modes appearing in the 
Lorentzian formulae of [14] are related to the Euclidean i? n -modes by r n = iR n . Moreover, 
for flat Euclidean backgrounds bounded by a 3-sphere, we choose the gauge-averaging term 
where $ is defined as (cf (5.1)) 

r)A °° °° 

d> = + ^VM, = J2 Rn(r)Q in \x) - r" 2 £ g n {r)Q^\x) . (3.1) 

n=l n=2 

Note that in (3.1) we have used the property of flat backgrounds ( 4 )V°Ao = and 
the relation s lJ = t~ 2 c 1 ^ between the contravariant 3- metric and the contravariant 
3-metric c lJ on a unit 3-sphere. Covariant differentiation on a unit S 13 , denoted by a 
vertical stroke, yields ' fc (x) = — Q^ n \x) [14]. With this choice of gauge-averaging 
functional, the corresponding differential operator acting on i? n -modes will turn out to 
be the one-dimensional Laplace operator for scalars if a = 1, as we would expect in the 
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light of the expansion (1.3a). Note also that (3.1) does not represent the Lorentz gauge- 
averaging functional (various alternative possibilities are studied in section 5). Thus, the 
part lE(g,R) of the Euclidean action quadratic in gauge modes is in our case [9] 



r / g n ■ n - 

+ 7r( \-rR n 

2a \ t 



dr (3.2) 



where we have inserted the flat-background hypothesis N = 1, a(r) = r. The physical 
degrees of freedom and the ghost field decouple from (3.2). Because we are here dealing 
with all degrees of freedom, we need further boundary conditions on the modes for Aq 
and the whole of A k . For example, we may set to zero on S 3 the whole of A k : f n (l) = 
<7n(l) = 0, Wi > 2. This, of course, implies the vanishing on S" 3 of the magnetic field 
B, whereas the converse does not hold, because B only depends on the / n -modes. As 
explained in [17], in this case the gauge-averaging term has to vanish as well on S" 3 by 
virtue of Becchi-Rouet-Stora-Tyutin (BRST) invariance. In light of (3.1), this implies that 
R n (l) = 0, Wi > 1. The ghost operator is then self-adjoint only if Dirichlet boundary 
conditions are imposed. Viceversa, if Neumann boundary conditions are chosen for the 
ghost field, remaining boundary conditions compatible with BRST invariance are / n (l) = 
and <7n(l) = 0, Vn > 2; R n (l) = 0, Vn > 1. This case is then called electric. 

We now integrate by parts in (3.2) and use the generalized magnetic or electric bound- 
ary conditions described above. Thus, defining Vn > 2 the second-order differential oper- 
ators 

An{T) = -J^-rTr + ^^ (3 - 3) 
13 
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A _ If <? 3 d \ (n 2 - 1) 



(3.4) 



we find Vn > 2 the fundamental result [9] 



r(n) — 1 



, (A n # n )dr + ^-/ T S R n (B n R n )dr 

* Jo \ n ~~ J-J ^ 



I/ 1 



+ 1- 



Tfif n i? n dr + / # n -Rn dr . 



(3.5) 



By virtue of gauge-invariance, we can perform the £(0) calculation setting a = 1, so that 
the contribution of rg n R n vanishes. However, our problem remains a coupled one, also 
with this choice. The eigenvalues for g n -modes and i? n -modes can be obtained in principle 

from the boundary conditions and the variational principle 5 (Je — A Je j = 0, where 

J E = \ J A a A^ Vdet g d 4 x. We here make the analytic continuation to the Euclidean time 
variable r = it in computing Ie (cf (3.2)) and Je, and we use (1.3)-(1.5) and the well- 
known properties of longitudinal vector harmonics and scalar harmonics [14] . This leads 
to the coupled system of two second-order ordinary differential equations for arbitrary a 
and Vn > 2 (the case n = 1 only involves the .Ri-mode, and should be treated separately, 
as in section 4) 



(n 2 - 1) 



9n 



I) 



(XT 



+ 1 



1 



OL 



t R n + R n — 



(n 2 - 1) 



rg n 



(3.6) 



— ( — Rn Rn J + o ~R n 



a 



r f j„ ( 1 . - -) + ^ = A„ ; '•/,>„ 
a / a 



(3.7) 
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Now we still choose a = 1, because it enables one to decouple much more easily the system 
(3.6)-(3.7). The boundary conditions are regularity at the origin 



g n (0) = R n (0) = foralln>2 



(3.8) 



and magnetic conditions on S 3 



g n (l) = R n (l) = for all n> 2 



(3.9) 



or electric conditions on S 3 



g n (l) = R n (l) = for all n > 2 . 



(3.10) 



In the a = 1 gauge we can express R n from (3.6) as 



n ~ t 1 i \ T 9n * ( o i\ 

[n z — 1) [n z — 1) 



'.. . 9n (n 2 -l) 

<?n H o 



(3.11) 



and its insertion into the corresponding form of (3.7) yields the fourth-order equation 



= 



~ (n z -l)J n (n 2 -l) 



W-(, 2 -i) 



+ 2r 



3A n r 2 
(n 2 - 1) 



9n + 



2r 2 



(n 2 - 1) 



n 



4 - \ n T< 



6-r 



(n 2 -iy n (n 2 -l) 



9™ 



(3.12) 



Moreover, studying first the magnetic case, the relations (3.9) and (3.11) lead to 



Sn(l) 0^(1) 



for all n > 2 



(3.13) 
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Of course, as shown by (3.6)-(3.7), the eigenvalues \ n have dimension (length) -2 . However, 
in (3.13) we have set a = 1 for simplicity, following (3.9)- (3. 10). Hence the physical 
dimension does not appear explicitly. For the solutions of the equations (3.11)-(3.12) 
subject to the boundary conditions (3.8)-(3.9), an existence and uniqueness theorem holds. 
Thus, denoting by k an integer > 0, in the light of the form of (3.12) we write its solution 



The insertion of (3.14) into (3.12)-(3.13), and the requirement that # n (l) = 0,Vn > 2, leads 
to a problem formulated in purely algebraic terms. One then finds that only half of the 
Q>n,k coefficients are non-vanishing and obey very involved recurrence relations, whereas 
the value of fx is obtained by solving a fourth-order algebraic equation. 
In fact, defining 



as [9] 



oo 




(3.14) 




(3.15) 



we find 




(3.16) 



F(m,n,n)a nim + 



2 - 



2(m + n) 2 -4(m + n) + 3 



) 



A n Q.n,m-2 — 



(3.17) 
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where m = 2,3, whereas, Vfc > 4, we have 

(2(/e + /i) 2 -4(fc + M)+3) 



= F(k,n,fi)a n , k + 



2 - 



(n 2 - 1) 



k-2 



(n 2 - 1) 



Q>n,k— 4 



(3.18) 



In (3.16)-(3.18), the value of /x can be obtained from the equation F(0,n, //) = 0, and 
bearing in mind (3.8), which implies that only a \x > 1 is an acceptable value, in the light 
of (3.11). In other words, we study the fourth-order algebraic equation [9] 



At 4 - (2n 2 - l)^ 2 + (n 2 - l) 2 = 



(3.19) 



This equation can be easily solved setting /j 2 = x and studying the corresponding second- 
order equation for x. One thus finds the four roots 



(3.20) 



,<?> = ^ - 1 - \ 



(3.21) 



(3.22) 



„L 2 > = -»? 



(3.23) 



Interestingly, both and ^ are > 1, Vn > 2. They yield the desired regular solution 



of the system (3.6)-(3.7). 
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4. Decoupled gauge mode 

In the light of (3.2), if we set a = 1 and integrate by parts using (3.8)-(3.10), we find that 
the decoupled gauge mode Ri(r) obeys the eigenvalue equation 

R 1 + -R 1 + fi 1 R 1 = (4.1) 
r 

which is solved by Ri(t) = Hit -1 J\(yfjl{ r), where Hi is a constant. Thus, in the 
magnetic case, we study the eigenvalue condition (setting a = 1 for simplicity) 

Ji (v 7 ^) " Ji (v 7 ^) = (4-2) 
whereas in the electric case the corresponding eigenvalue condition is 

Ji(v^T) = . (4.3) 

If (4.2) holds, we have to perform a ((0) calculation with just one perturbative mode, 
subject to a complicated eigenvalue condition involving a linear combination of J\ and J\ . 
Of course, a regularization is still needed because there are infinitely many solutions X n 
of (4.2). For this purpose, it is convenient to use the technique described and applied in 
[3,9,18]. The basic idea is as follows. 
Given the C-fimction at large x 

oo _ 

C(*,x 2 )^(\ + * 2 )" S (4.4) 

n=l 
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one has in four-dimensions 

r(3)C(3,x 2 )= / tV l2t G(t) ( it~VB ? r(i + ;)x-'- 2 (4.5) 

J ° q =o v z/ 

where we have used the asymptotic expansion of the heat kernel G(t) for t — > + , written 



as 



Such an asymptotic expansion does actually exist in the case of Laplace operators subject 
to Dirichlet, Neumann or Robin boundary conditions [9,13,19,20]. On the other hand, 
defining (cf (4.2)) 

F 1 (z) = J 1 (z)-zJ l (z) (4.7) 

one also has the identity 

T(3)C(3,x 2 ) = -iV 1 (-^^ 3 log^i*)- 1 ^!*)) (4.8) 

where N\ = 1 is the degeneracy of the problem. Thus, the comparison of (4.5) and (4.8) 
can yield the coefficients B q and in particular £(0) = B4, provided we carefully perform an 
uniform asymptotic expansion of F\(ix). Now, making the analytic continuation x — > ix 



and then defining a = yl + x 2 , one obtains the following asymptotic expansions which 
are uniformly valid in the order as | x |— > 00: 

Mix) ~ ^=La-h«e- ^(^E^l, a(x)) (4.9) 

V 27T 
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J\(ix) ~ -^^eV log(1+;) E 2 (l,5(a;)) (4.10) 
V2vr 

where Ei(l,5(x)) ~ ££L S 2 (l,5(x)) ~ ££W£) [3,9,18]. Using (4.7) and 

(4.9)-(4.10), and defining C = — log(\/27r), one thus finds the uniform asymptotic expan- 
sion 



oo I 

log^ix)- 1 ^^)) ~C-log(l + 5) + -log(5)+5 + ^^6 Zr 5- z - 2r (4.11) 



i=l r=0 



where the double sum on the right-hand side of (4.11) is obtained by expanding in in- 
verse powers of a the log^S*- — £ 2 j. In the light of (4.8) and (4.11), we conclude that 

r(3)C(3, x 2 ) ~ (7i + <r 2 ] , where 



1 d 
o~\ ~ l — 

2x ax 



log(l + a) + - log(a) + a 



(4.12) 



0"2 



2x dx 



3 oo / 



EE" 

1=1 r=0 



-Z-2r 



(4.13) 



It can be easily checked that the asymptotic expansion of o<i in (4.13) does not contribute 
to C(0), whereas one finds 



3 _ x~ 6 ^ 

o-i ~ -x" 5 - x~ 6 + —— + 2_ J u k x~ k 



8 



(4.14) 



k=7 



which implies (cf (4.5)) 



(4.15) 
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By contrast, when the eigenvalue condition (4.3) holds for the i?i-mode, the square bracket 
on the right-hand side of (4.12) contains — |log(5) rather than | log(5) (cf (4.9)). This 
leads to 

C*,(0) = i( -!)=-! (4.16) 

since, again, the asymptotic expansion of the corresponding o~i does not contribute to 
C(0). Note also that our value (4.16) agrees with the result found in [21] in the case of 
one perturbative mode subject to Dirichlet boundary conditions. As explained in [22], the 
perturbative calculations of [21] are correct, whereas the ((0) values with finitely many 
perturbative modes of appendix A of [23] are incorrect. 

5. Ghost-field contribution 

In section 3, we have chosen the gauge-averaging term where the gauge-averaging 
functional $(^4) can be written in the form (cf [24]) 

$(A) = ee WVA, - K\A = ^ + (3) V*A, (5.1) 

OT 

where K\ = ^ is the trace of the extrinsic-curvature tensor of the boundary. This choice 
of &(A) leads to the familiar one-dimensional Laplace operator acting on the i? n -modes, 
which simplifies the £(0) calculation for the coupled gauge modes and for the i?i-mode, 
as shown in sections 3 and 4. However, since $i(A) is not the Lorentz gauge-averaging 
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functional, the corresponding ghost action does not involve the familiar Laplace operator. 
This is proved (cf [15]) by studying the gauge transformation 

^ = ^ + ^ = ^ + 9,6 (5.2) 

where the scalar e is expanded on a family of 3-spheres centred on the origin as 

oo 

e(x,r) = J2^n(r)Q in \x) . (5.3) 

n=l 

One thus finds 

oo 
n=l 

This implies that the eigenfunctions of the ghost operator are of the kind [16] 

e n {r) = V^J^ l {^Er) . (5.5) 

More precisely, since the electromagnetic field is bosonic, the corresponding ghost field 
is fermionic [15]. Its contribution to the full ((0) is thus obtained changing sign to the 
scalar-eigenfunctions contribution of (5.5), and then multiplying the resulting number by 
two, since the ghost field is complex. We now have to perform a ((0) calculation which 
involves Bessel functions of non-integer order, generalizing the technique described in sec- 
tion 4. Here we show that, although eigenvalues and eigenfunctions are different, the ((0) 
calculation originating from (5.5) is closely related to a standard ((0) calculation involving 
Bessel functions of integer order. For this purpose, we study the simplest case, i.e. when 
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d 2 (n 2 - 1) 



dr 2 



e n(T) 



(5.4) 
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the ghost field obeys homogeneous Dirichlet conditions on S" 3 . This leads to the eigenvalue 
condition 

J^^(VEa) = for all n > 1 . (5.6) 
Following [9], our section 4 and (5.6), it is now useful to define Vn > 1 and at large x 



v = +\ n A — 



(5.7) 



a v {x) = a/z/ 2 + x 2 = ^2 - | + 



a n (x) = \/n 2 + x 2 



(5.8) 



(5.9) 



Since the generalization of the technique of section 4 to infinitely many perturbative modes 
for the ghost involves defining a v (x), whereas we are only able to perform exact calculations 
using a n (x), it is also useful to evaluate the ratio 



Oiy(x) 

a n (x) 



Pn(x) 



1 - | (n 2 + x 2 ) _1 - A (n 2 + x 2 ) " 2 + O ( (n 2 + x 2 " " 



(5.10) 



The asymptotic expansion (5.10) is very useful in that it is uniform in n, i.e. it holds Vn > 1, 
at large x. A careful study of section 7.3 of [9] shows that, if the eigenvalue condition (5.6) 
holds (whose eigenvalues are positive Vn > 1), (4.8) and (4.11) are generalized as 



T(3)C(3,x 2 ) 



(5.11) 
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where 



E" 2 

n=0 



ux~ 6 + p 2 x~ 6 a~ 1 + V —x 4 a„ 3 + \v 2 x 2 a~ 5 - ^- + \a v 5 

2 8 2 8 



(5.12) 



OO I 



^2 ~ " ai - ( r + 2 ) ( r + 2 + *) ( r + 2 + 2 ) ^ n2p2ra - 

1 = 1 r I) 1 " " ' n=0 



l-2r-6 



(5.13) 



In these formulae, obtained using uniform asymptotic expansions of Bessel functions of 
non-integer order, n 2 is the degeneracy resulting from the scalar harmonics appearing in 
the expansion (5.3), v is the order of the Bessel functions defined in (5.7), and a u has been 

defined in (5.8). We now re-express v 2 as (n 2 — |j, and a v (x) ~ a n (x)p n (x) as in (5.10). 

Moreover, we use the contour formula [3,9,18] 



E 

n=0 



. r(* + Qr(?-j) 



2r 



0+f) 



for all k = 1,2,3, 



(5.14) 



We then point out that the asymptotic expansion (5.12) can be cast in the form 



Ol 



- x-*l£ + x-H£ + x" 4 4 3 ) + x~ 2 l£ - l£> + /<?] (5.15) 



where (see appendix A) 



OO OO 



n . n ."f V 7 



n=0 n=0 



(5.16) 
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1 °° Q °° inr 00 

00 2 ^ n 16 *» 256 ^ 

n=0 n=0 n=0 

3 ^ — \ 9 — s 27 v 9 — ^ 405 ^ — \ 2 —7 

> na^ > n a„ > n a ' 

8^ n 64 n 1024 ^ 

n=0 n=0 n=0 

1 00 „ 00 

+ ^E- 4 «n 3 0(a n - 6 ) -^n 2 a n - 3 0(a n - 6 ) 

n=0 n=0 

JCJ) „ I y n 4 a -5 45 y 4 7 JM5 y 4^-9 
00 g n n 1024 ^-^ 

n=0 n=0 n=0 

9 v -v 2 —5 135 v -v 2 —7 2835 v -v 2 — q 

> n a„ > n a„ — > n a„ 

32 ^ n 256 ^ n 4096 ^ 

n=0 n=0 n=0 

„ 00 Q 00 

+ i E - 4 «n 5 0(« n - 6 ) n 2 a n - 5 0(a n " 6 ) 

n=0 n=0 

I<?^f> 2 a- fl 

T1=0 

-.OO „ OO ^ OO 

- ^ E - 2 «n 6 + 1 E - 2 «n 8 + ^ E - 2 «n 6 o (V ) 

n=0 n=0 n=0 

n=0 

3 00 45 00 945 °° 

~ - n 2 a~ 5 H n 2 «~ 7 H n 2 «~ 9 

8^ n 64 ^ n 1024 ^ 

n=0 n=0 n=0 

+ ^E^ 5 0(«; 6 ) • 

n=0 
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It is therefore clear, using (5.14), that the £(0) value resulting from o\ and cr 2 is given 

by = ~ 2 ( — yIo) ( wn ich coincides with the C(0) value corresponding to the Lorentz 

gauge-averaging functional) plus additional terms owed to the second sum in (5.16), the 
third and fifth sum in (5.17)-(5.18), denoted by Ti, T 2 , T 3 , T 4 , the third sum in (5.20), 
denoted by T 5 , and finally (5.13). Note that defined in (5.19) does not contribute to 
the additional terms. The detailed calculation yields 



n=0 



^ 2 --S-- 4 f;- 2 ^ 5 = -^- 6 (5-22) 

n=0 



n=0 



n=0 

T 5 = ^-Yn 2 a- 9 = —x- e . (5.25) 
5 1024 ^ n 128 K ' 

We now focus on (5.13) and (5.16), and we first study the asymptotic expansion (5.13), 
since (5.16) gives rise to severe technical difficulties (see below). For this purpose, we 
remark that, studying for all integer values I e [1, oo[, r e [1, I] the function (see appendix 
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A) 



n=0 



l-2r-6 



J2 n 2 (n 2 - l) r a~ 1 -^ 1 + A lr a~ 2 + B w ctf + o(o- 6 ) 



n=0 



r(l) , r(2) , A3) A4) 
A lr "r J Zr "r J Zr J Zr 



one finds 



oo 

i 1} (-)-E- 2 ( 



n=0 



/< 2 >(*)~^£n 2 (n 2 -?) 



a. 



n=0 



/« 3 »(*)~B !r 5> 2 („ 2 -2) 



a. 



-Z-2r- 8 



-Z-2r-10 



n=0 



00 „ r 

i 4) (-)-E- 2 (- 2 "l) «n^- 6 0(a n - 6 ) 



n=0 



(5.26) 



(5.27) 



(5.28) 



(5.29) 



(5.30) 



where A\ r and B\ r are coefficients which only depend on I and r. The case r = is easier. 
Using (5.13)-(5.14), r = leads to a contribution to £(0) related to 



T 6 = — aio 



15 21 



<(§)r(3) 



8 8 2r ( f ) 



x 



--aiox 



(5.31) 



where -g- is the coefficient of a n 2 in the asymptotic expansion of p n 7 (x) (see appendix A). 
If the integer r is > 1, one has to study (5.26)-(5.30), where 



("•-§)' 



2r ^ 2r-2 , , 2 ( 3 

n — -rn + ... + rn 
4 



( 3 \ r_1 / 3 V 



(5.32) 
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Inserting (5.32) into (5.27)-(5.30), and using (5.14), a lengthy calculation yields a contri- 
bution to C(0) related to (see appendix A) 



T 7 



3 9 
4 an " g an 



x 6 = —-anx 6 



(5.33) 



Note that the two terms on the r.h.s. of (5.33) are due to the asymptotic expansions of 

l\l\x) and l\l\x) respectively, whereas (5.29)-(5.30) do not affect the £(0) value, since 
they do not involve x~ 6 . It now remains to evaluate the contribution of (5.16). Indeed, 
defining 



OO , x 

n=0 x ' 



(5.34) 



we point out that multiplying and dividing the round bracket by {i>-\-a u j , and then adding 
and subtracting a u in the numerator of the corresponding expression, one finds by virtue 
of (5.8) the useful identity 



n 



n=0 



- /(I) + 7(2) 

— J oo ' OO 



(5.35) 



Moreover, (5.10) and (5.14) show that the contribution to ((0) owed to is related to 



128 



-x 



n=0 



128 



(5.36) 



A further contribution is owed to 



T ^Wo x 
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originating from Y^=o 71,3 in (5-16). However, we do not yet know how to deal properly 
with the divergent sum 

„=o [y + ol v \ 

We should now add up the numerical coefficients appearing in (5.21)-(5.25), (5.31), 
(5.33), (5.36)-(5.37), divide them by two, and finally multiply by —2 since the ghost is 
fermionic and complex. This leads to the following partial contribution to the £(0) value 
for the ghost field: 

>m/^ 19 3/ \ 9 1 1 11 

C ^ (0) = 90 - 512 + 8 I" 10 + H + T28-T20 = 360 + 512 (5 " 39) 

where ^ is added for the reasons described following (5.20), and we have used the values 
aio = g, an = — ^ appearing in equation (26) of [18]. 

By contrast, if the Lorentz gauge-averaging functional is chosen, one finds 

§(A) = $ 2 (A) ee WVi, = ^ + (4) VM, (5.40) 

which implies 

oo 

8($ 2 (A)) ee 0> 2 (A) - ^A) = Q in) (x) 

n=l 

where we have used the property ^V^e = ^V^e = = die, Vz = 1,2,3. Thus, as we 
anticipated, the familiar one-dimensional Laplace operator appears in the ghost action, so 
that the ghost contributions to the full £(0) value are more easily computed as —2 (— -jgo) 
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dr 2 t dr 



3 d (n 2 
+ — 



1) 



e n (r) (5.41) 
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and —2 (j§q) in the Dirichlet and Neumann cases, respectively. However, if $2(^4) is 
chosen as gauge-averaging functional, the form of the action quadratic in the gauge modes 
becomes Vn > 2 



2 Jo (n 2 -l) 



d 2 g n 1 alg n | (n 2 - 1) 
dr 2 r dr ar 2 ® n 



dr 



+ 



2io 



r 3 R n 



d 2 R n 3 dR r 



dr 2 



+ 



r dr 



+ n 2 -l + 



CK/ 7" 



dr 



1 



a 



+ 1 / r 9nR n dr+ 1 / g n R n dr 



a 



rg n R r , 



i 1 
+ — 

o 



2a 



T 3 R n Rr 



(5.42) 



Thus, the second-order differential operator acting on i? n -modes is no longer the one- 
dimensional Laplace operator for scalars, and the calculation becomes more involved. For 
example, if we set a = 1, the contribution of -Ri(r) to ((0) involves a Bessel function of 
order Moreover, a no n- vanishing boundary term 1^ = ^R n (a)R n (a) = — ^i? 2 (a) 
survives in the action, if the whole functional $2(^4) is required to vanish on the boundary 
in the magnetic case (cf [17]). Thus, one has to add to the action a boundary term equal 

to — , if the whole of $2(^4) is set to zero on S" 3 . 

Of course, since the theory is gauge-invariant, infinitely many other choices for <fr(A) 
(but not all choices) are still possible. A very relevant class of choices can be cast in the 
form 

$W (4) = (4) + bK\A (5.43) 

where b is a real number. With our parametrization, b = — 1 leads to $i(A), and 6 = 
leads to $2(^4)- Note that, even if we set a = 1, it does not seem possible to decouple 
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gauge modes using |^ and obtain a well-defined ghost action, since the decoupling of g r 
and R n , Vn > 2, is then obtained setting 



*(A) = $ 3 (A) = J] a 

n=2 \ 



Q {n \x) . (5.44) 



However, the ghost action should be derived by functionally differentiating the infinite sum 
of square roots on the right-hand side of (5.44) as in (5.4) and (5.41), and this does not lead 
to a linear, second-order differential operator. This is why we believe that the coupling of 
gauge modes is an intrinsic property of problems with boundaries, as well as the choice 
of gauge-averaging functionals of the form (5.43), which all reduce to the Lorentz choice 
in the absence of boundaries. Note that the work in this section supersedes earlier work 
appearing in section 6.5 of [9], where the ghost-field operator (cf (5.4)) was not derived. 

In light of (2.6), (4.15) and (5.39), the full ((0) value for vacuum Euclidean Maxwell 
theory in the case of magnetic boundary conditions on S" 3 takes the form 



C(o) = C ( b PDF) (0) + CrM + Cgm(0) + ( gh (0) = ~ + ii + Cgm(0) + C^(o) (5.45; 



where Cgm(0) and Cgl^(O) are the as unknown contributions to ((0) arising from 
coupled gauge modes (section 3) and from (5.38) respectively. We have been unable to 
evaluate Cgm(0) since we do not know explicitly the uniform asymptotic expansion as 
A n — > oo of the power series in (3.14), which is not (obviously) related to well-known 

special functions (see appendix B). Moreover, the regularized contribution (0) of (5.38) 
to C(0) involves v = +\Jn 2 — |, which is a source of complication. However, it should be 
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emphasized that all divergences are only fictitious, since the starting point for the derivation 
of (5.12) is the identity [9] 



2x dx J V v + a v 



= [V + OL v 



—•x 9 _a 3 2 _ ^ 

a., va„ v a„ 

8^8^ 



(5.46) 



This proves that by summing over all integer values of n from to oo one gets a convergent 
series. 

In this section we have not studied the case of Neumann boundary conditions for the 
ghost field, i.e. the electric case. This complicated calculation may be, by itself, the object 
of another paper. However, interestingly, in light of (2.6), (2.25) and (4.15)-(4.16) one 
finds 



C { b DF) (0) + C* (o) = C1T ' ; (o) + ( Rl (o) = - 



(PDF). 



61 

90 



(5.47) 



In other words, if the gauge-averaging functional of (5.1) is chosen, physical degrees of 
freedom and decoupled gauge mode give the same partial contribution to the full C(0), i-e. 



— fjj, both in the magnetic and in the electric case. 



6. Concluding remarks and open problems 



One-loop quantum cosmology may add further evidence in favour of different approaches 
to quantizing gauge theories being inequivalent [2-9,25-34]. Studying flat Euclidean back- 
grounds bounded by a 3-sphere, for vacuum Maxwell theory the PDF method yields 
C(0) = —j^q and C(0) = j§q in the magnetic and electric cases respectively [6,9,14], 

32 



Gauge-averaging junctionals for Euclidean Maxwell theory ... 
whereas the indirect method (by this we mean that one-loop amplitudes are expressed 
using the boundary-counterterms technique and evaluating the various coefficients in a 
covariant way as in [2,17]) was found to yield £(0) = — f§ in both cases in [2]. For N = 1 
supergravity, the PDF method yields partial cancellations between spin 2 and spin | [7-9] , 
whereas the indirect method yields a one-loop amplitude which is even more divergent than 
in the pure-gravity case [2]. Finally, for pure gravity, the PDF method yields £(0) = — ^ 
in the Dirichlet case, whereas the indirect method yields £(0) = — ^ [2,9]. Moreover, 
within the PDF method, it is possible to set to zero on S 3 the linearized magnetic cur- 
vature. This yields a well-defined one-loop calculation, and the corresponding £(0) value 
is ^ [9]. By contrast, using the Faddeev-Popov formula, magnetic boundary conditions 
for pure gravity are ruled out [2]. Interestingly, recent work in [35] seems to add evidence 
in favour of direct £(0) calculations being correct. The authors of [35] have shown that 
different formulae for £(0) previously obtained in [8] for Majorana and Dirac fermions on 
the part of a de Sitter sphere bounded by a 3-sphere, with local and spectral boundary 
conditions, have the same limiting value in the case of a full sphere. This value coincides 
with the covariant one obtained by the method in [2,13,17]. Moreover, all these expressions 
are found to give the same results in the case of a hemisphere [35]. The authors of [35] 
have also suggested that 3 + 1 splits of the kind considered in many papers, including Eqs. 
(1.3a)-(1.3b) of our paper, might be the reason of the discrepancies for higher-spin £(0) 
values found using covariant and non-covariant methods. However, there is not yet a proof 
of this statement, and in the case of real scalar fields subject to Neumann (or Dirichlet) 
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conditions on a 3-sphere boundary, £(0) values obtained from various methods coincide, 
although the boundary 3-geometry is the same as in higher-spin calculations. 

It is therefore necessary to get a better understanding of the manifestly gauge-invariant 
formulae for one-loop amplitudes used so far in the literature, by performing a mode-by- 
mode analysis of the eigenvalue equations, rather than relying on general formulae which 
contain no explicit information about degeneracies and eigenvalue conditions. This detailed 
analysis has been attempted here in the simplest case, i.e. vacuum Maxwell theory at one- 
loop about a flat Euclidean background bounded by a 3-sphere. Our results are here 
summarized for the sake of clarity: 

(1) In the light of (2.6), (2.25) and (4.15)-(4.16) the physical degrees of freedom, and 
the decoupled gauge mode, give a contribution to the full £(0) equal to — |g both in the 
magnetic and in the electric case (this important property had not been realized in section 
6.5 of [9]), if the gauge-averaging functional <&\{A) of (5.1) is chosen. Since in [2] it was 
found that the full C(0) values for spin 1 are equal in the magnetic and electric cases, it 
appears relevant that also our partial contributions to the full £(0) coincide in these two 
cases for the spin-1 problem about flat Euclidean backgrounds. 

(2) Remaining gauge modes g n and R n always obey a coupled system of linear, second- 
order ordinary differential equations, Vn > 2. The solution of such a system corresponding 
to &i(A) has been given in section 3 and appendix B. 

(3) If &i(A) is chosen, the ghost eigenfunctions involve Bessel functions of non-integer 
order. The corresponding contribution to the full £(0) can be obtained using the method 
of section 5. Such a technical point appears interesting, since to our knowledge no previous 
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mode-by-mode analysis for the ghost is appearing in the literature in the case of Bessel 
functions of non-integer order. 

It now remains to evaluate the contribution to the full £(0) of the divergent sum in 
(5.38), and the uniform asymptotic expansion of g n - and i? n -modes as X n — > oo at the end 
of section 3. Unfortunately, the generalization of the method described in [5-8] is highly 
non-trivial. By contrast, a simpler form of the ghost eigenfunctions is obtained using the 
Lorentz gauge-averaging functional <3>2 (A). However, this leads to a further complication of 
the calculations involving gauge modes, since the decoupled mode Ri(r) involves a Bessel 
function of order y/l3 (this implies a contribution to ((0) proportional to VT3, which we 
find very puzzling), and coupled gauge modes require the addition to the action, in the 

magnetic case, of a boundary term equal to ^ ^^=2 -^n( a )- 

Thus, although some evidence exists that different £(0) values for gauge fields in the 
presence of boundaries are due to inequivalent quantization techniques [9] , the most im- 
portant check, i.e. the mode-by-mode analysis of eigenvalue equations for gauge modes 
and ghost fields, remains a very difficult problem. We hope that our paper, through its 
detailed (although incomplete) analysis, may contribute to shed new light on this long- 
standing problem in quantum field theory. 
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Appendix A 

The derivation of (5.17)-(5.20), (5.31), (5.33) and (5.36) has been obtained using the 
following asymptotic expansions: 




(Al) 
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Pn\x) ~ 1 + f < 2 + ^«; 4 + 0(V) (AT) 

p- 8 (s) ~ 1 + 3a~ 2 + y«„- 4 + 0(a" 6 ) (A8) 

p; 9 W-l + f«; 2 + ^|«; 4 + 0(a; 6 ) • (A9) 

Note that these expansions are valid uniformly in the integer n, Vn > 1, as | x | — > oo. 
They are obtained using repeatedly (5.10) and the well-known expansion of (1 + Y) _1 as 
Y ->0. 



Appendix B 



Following section 3, coupled gauge modes can be written as 



oo 

#(r) = E a S(^Ml 3 ' ) )^ (5-1) 



OO 

R^\r) = ^ k {n,k,X^)r k ^- 1 . (5.2) 



The label j is introduced because, for each integer value of n > 2, there is a whole family 
|An^| of eigenvalues labelled by the integer j, say. They are solutions of the equation 

9n( a ) = 0, and their degeneracy dj(n) = n 2 , Vj > 1 and Wn > 2 [14]. Now, defining Wk > 2 
and Wn > 2 

f2(A; + ^) 2 - 4(fc + /x) + 3) 

n, //) = 2 - A L (jB . 3 ) 

— 1J 
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one finds Vj > 1 

,00 



a 



n , 2 G(2,n,fi) 



Aj?) (5.4) 



Ai j) ) 2 

F(fc, n, ^ogi + G(fc, n, /*) A^agU - ^ifZ^U = Vfc > 4 (5.5) 

iVS (*-6) 



y n,0 



(n 2 - 1) 



a "i = 6 "i = Vfc = (2m + 1) m = 0,l,2,... . (fl.8) 

Moreover, setting the 3-sphere radius a to 1 for simplicity, magnetic boundary conditions 
(i.e. g n (l) = R n (l) = 0) lead to 

oo 
A;=0 

v^°° 1,3 U) J.2 (j) 

A£> =fn 2 - l) - ^ - 2) - ^ fc= ° ^ - (3a* - 1) ^ fc= ° • (5-10) 

Since, Vn > 2, there are two values of [i > 1, a further label is necessary to characterize 
completely the coupled gauge modes as follows: 

9$ («> Agi, r) and flft (n, r) if ^ = 
9 { 2,l (n, Agi, r) and (n, A^, r) if p = ^ 
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(see (3.20)-(3.21)). 

Note that it is extremely difficult (if not impossible) to find the eigenvalues An ^ by 
analytic or numerical methods, since (B.4)-(B.5) imply that a function H exists such that 

-^ = H(k,n,»)(Wy (B.ll) 

°n,0 

for all even values of k > 2, and Vn > 2. Thus, when (B.ll) is inserted into (B.9)-(B.10), 
it is not clear how to find an explicit solution for and a^\. 
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